FAST SOLITON SCATTERING BY ATTRACTIVE DELTA 

IMPURITIES 
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^P ' Abstract. Wc study the Gross-Pitacvskii equation with an attractive delta func- 

tion potential and show that in the high velocity limit an incident soliton is split 
into reflected and transmitted soliton components plus a small amount of disper- 
"gj ■ sion. We give explicit analytic formulas for the reflected and transmitted portions, 

["t I , while the remainder takes the form of an error. Although the existence of a bound 

f^ ' state for this potential introduces difficulties not present in the case of a repulsive 

potential, we show that the proportion of the soliton which is trapped at the origin 
vanishes in the limit. 
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1. Introduction 
The nonlinear Schrodinger equation (NLS) or Gross-Pitaevskii equation (GP) 



CsJ ■ (1-1) idtu + Idlu + \u\'^u = , 

> ' 

lO ! where u = u{x, t) and x G M, possesses a family of soliton solutions 

m ' u{x, t) = e'^e™"e-^'*"'Asech(A(x - xq - vt)) 

O ; 

— u ', parameterized by the constants of phase 7 G M, velocity f G M, initial position xq G 

O I and scale A > 0. Given that these solutions are exponentially localized, they very 

(^ [ nearly solve the perturbed equation 

(1.2) idtu + ^dlu - q6o{x) + \u\\ = , 

when the center of the soliton |xo + t't| ^ 1. In fact, if we consider initial data 

u{x, 0) = e*'"'sech(x - Xq) 

for xq ^ — 1, then we expect the solution to essentially remain the rightward prop- 
agating soliton e*^''e~2™ *sech(x — Xq — vt) until time t ~ |xo|/f at which point a 
substantial amount of mass "sees" the delta potential. It is of interest to examine the 
subsequent behavior of the solution, as this arises as a model problem in nonlinear 
optics and condensed matter physics (see Cao-Malomed [1] and Goodman-Holmes- 
Weinstein [2]). In the case |g| <^ 1, Holmer-Zworski [5] find that the soliton remains 
intact and the evolution of the center of the soliton approximately obeys Hamilton's 
equations of motion for a suitable effective Hamiltonian. This result applies to both 
the repulsive (g > 0) case and the attractive (g < 0) case, and identifies the |g| <^ 1 
setting as a semi-classical regime. 
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On the other hand, quantum effects dominate for high velocities |i;| ^ 1. In Holmer- 
Marzuola-Zworski [3] [4] , the case of g > and t; ^ 1 is studied (most interesting is the 
regime g ~ w), and it is proved that the incoming sohton is spht into a transmitted 
component and a reflected component. The transmitted component continues to 
propagate to the right at velocity v and the reflected component propagates back to 
the left at velocity —v, see Fig. 1. The transmitted mass and reflected mass are 
determined as well as the detailed asymptotic form of the transmitted and reflected 
waves. The rigorous analysis in [3] is rooted in the heuristic that at high velocities, 
the time of interaction of the solution with the delta potential is short, and thus 
the solution is well-approximated in L^ by the solution to the corresponding linear 
problem 



;i.3) 



idtu + \d1u — q5o{x)u = . 



This heuristic is typically valid provided the problem is L^ subcritical with respect 
to scaling. In this case, it is shown to hold using Strichartz estimates for solutions 
to this linear problem and its inhomogeneous counterpart, with bounds independent 
of q. The Strichartz estimates are also used in a perturbative analysis comparing the 
incoming solution (pre-interaction) and outgoing solution (post-interaction) with the 
solution to the free NLS equation (1.1). One then proceeds with an analysis of the 
linear problem to understand the interaction. Let Hg = —^d^ + qSo{x) and consider 
a general plane wave solution to {Hg — |A^)w = 0, 



w{x) 



A^e 



~i\x 



B^e'^"" for X > 



iXx 



^_g-.Ax ^ D e*Ax for a; < 



The matrix 



SW: 






A_ 
B_ 



sending incoming (+) coefficients to outgoing (— ) coefficients is called the scattering 
matrix, and in this case it can be easily computed as 



S{X) 



tg{X) rg{X) 

MX) tg{\)^ 



where tg{X) and ^^(A) are the transmission and reflection coefficients 



t,{X) 



iX 



and rg{\) 



i\ — q 



iX — q 

We have that at high velocities and for Xi <^ —1, 

(1.4) 

e-'*^«[e'""sech(x - xi)] ^ t(w)e-**^o[e"''sech(x - xi)] + r(t;)e-**^°[e-""sech(x + xi)]. 
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From this we can infer that the transmitted mass 



^'W = "ww = ^""'""' 



matches the quantum transmission rate at velocity v, i.e. the square of the transmis- 
sion coefficient 

»,2 



T,{v) ^ Uv)\^ 



y2 _|_ q2 



This is confirmed by a numerical analysis of this problem in Holmer-Marzuola-Zworski 
[4], where it is reported that for q/v fixed, 

Tq(v) = ^5 5" + 0(v ^), as V ^ +00 . 

yZ _|_ qZ 

Further, (1.4) gives approximately the form of the solution just after the interaction, 
and one can then model the post-interaction evolution by the free nonlinear equation 
(1.1) and apply the inverse scattering method to yield a detailed asymptotic. The 
results of [3] are valid up to time logt', at which point the errors accumulated in the 
perturbative analysis become large. 

When g < 0, the nonlinear equation (1.2) has a one-parameter family of bound 
state solutions 

(1.5) u{x, t) = e**^'/2^sech(A|x| + tanh-^(|g|/A)), < |g| < A 

The numerical simulations in [4] show that at high velocities, the incoming soliton is 
still split into a rightward propagating transmitted component and a leftward prop- 
agating refiected component, although in addition some mass is left behind at the 
origin ultimately resolving to a bound state of the form (1.5). However, the amount 
of mass trapped at the origin diminishes exponentially as v ^ +cxd and the observed 
mass of the transmitted and refiected waves is consistent with the assumption that 
the outgoing solution is still initially well-modelled by (1.4). 

In this paper, we undertake a rigorous analysis of the g < and v large case. This 
analysis is complicated by the presence an eigenstate solution u{x,t) = ea**'' g-kll^l 
to the linear problem (1.3). Therefore, the Strichartz estimates, which involve global 
time integration, cannot be valid for general solutions to (1.3). However, they can 
be shown to hold for the dispersive component of the solution e~**^''((/) — P4>), where 
P is the orthogonal projection onto the eigendirection e~'^"^'. In the pre-interaction, 
interaction, and post-interaction perturbative analyses, this eigenstate must be sepa- 
rately analyzed. This introduces the most difficulty in the post-interaction analysis, 
although (as explained in more detail below), we are able to obtain suitable estimates 
by introducing a more refined decomposition of the outgoing waves and invoking some 
nonlinear energy estimates. We thus obtain the following: 
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Figure 1. Numerical simulation of the case q = —3, f = 3, xq = —10, 
at times t = 0.0,2.7,3.3,4.0. Each frame is a plot of amplitude |m| 
versus x. 

Theorem 1. Fix < e <^ 1. Ifu{x,t) is the solution of (1.2) with initial condition 
u{x, 0) = e*^^sech(x — Xq) and Xq < —v^, then for \q\ > 1 and 
(1.6) V > C(log \q\y/' + C|g|^(i+2.) ^ ^^^^^^^^i 

we have 



;i-7) 



\u{x,t)\'^dx 



2 .lx>o ■ ■ ■ ■ v- + q^ 



p. _|_ ^2 

uniformly for "post-interaction" times 



+ o{\q\h-'^(^-^'^) + 0{t 



-(l-2£)^ 



Fol 



+ V 



-1+s 



< t < erlogt;. 



Here the constant C and the constants in O are independent of q,v,e, while Ce,n is a 
constant depending on e and n which goes to infinity as e ^ Q or n ^ oo. 

The proof is outlined in §3. It is decomposed into estimates for the pre-interaction 
phase (Phase 1), interaction phase (Phase 2), and post-interaction phase (Phase 3). 
The details of the estimates for each of the phases are then given in §4 (Phase 1), §5 
(Phase 2), and §6 (Phase 3). 
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The assumption that v > |g|T4(^+2^-' is new to our g < analysis; no assumption of 
this strength was required in the g > case treated in [3]. It is needed in order to 
iterate over unit-sized time intervals in the post-interaction phase. The perturbative 
equation in that analysis has a forcing term whose size can be at most comparable 
to the size of the initial error. The condition that emerges is | gp/^ (error) ^ < c. 
Since the error bestowed upon us from the interaction phase analysis is \q\3v~6^^~^\ 
the condition | g|^'^ (error )^ < c equates to t; > |g| 14(^+2'^). Provided this condition is 
satisfied, we can interate over ~ e log v unit-sized time intervals, with the error bound 
doubling over each interval, and incur a loss of size v^. This enables us to reach time 
elogw. 

The assumption v > |g| 14(1+2'^) is not a serious limitation, however, since the most 
interesting phenomenon! (even splitting or near even splitting) occurs for |g| ~ v. 
Furthermore, if the analysis is only carried through the interaction phase (ending at 
time \xo\/v + v~^~^'^) and no further, then only the assumption v > |g|2(^+^) is needed. 
We believe that if our post-interaction arguments are amplified with a series of tech- 
nical refinements, we could relax the restriction needed there from v > |g|T4(i+2^) to 
V > |g|2(^+'^). On the other hand, the condition v > |g|2(^+^) shows up in a more 
serious way in the interaction analysis, and to relax this restriction even further (if it 
is possible) would require a more significant new idea. 

The condition v > |g| z^^^^-* comes about as a result of applying Strichartz estimates 
to the flow e~**^«0 rather than just to the dispersive part e~^^^''{(j) ~ P'P)-, and the 
additional error found in the g < case compared to the g > case arises in the 
same way. As discussed in Theorem 3 below, in the case g > we have 0{v~^^) 
in place of 0{\q\'^v~~>^) + 0{v~^'^) for g < (which in the crucial regime \q\ ~ v 
becomes 0{v~i^)). However, the numerical study conducted in [4] (see equation 
(2.4), Table 2, and Fig. 5 in that paper) suggests that the trapping at the origin 
should be exponentially small instead, indicating that this is probably only an artifact 
of our method of proof. 

The proof of Theorem 1 is based entirely upon estimates for the perturbed and free 
linear propagators, and some nonlinear conservation laws (energy and mass); there is 
no use of the inverse scattering theory. However, as in [3] , we can combine the inverse 
scattering theory with the proof of Theorem 1 to obtain a strengthened result giving 
more information about the behavior of the outgoing waves. This result we state as: 



Theorem 2. Under the hypothesis of Theorem 1 and for 



h 1 < t < elogf. 
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we have 
(1.8) 
u{x,t) = (/)o(|tq(tO|)e5^l^''(^)l'V'^''s*'('')e^^^e-**'''f;(t;)sech(f;(i;)(x - Xq - tv)) 

where 

(1.9) f,iv) = [2\t,iv)\ - 1] + , R,iv) = [2\r,iv)\ - 1]+ , 

0o(«) = y^ i-g^i+^^^KvJc+(2«-i)^'^ 

VK/ien2|tg(i;)| = 1 or2\rq{v)\ = 1 the first error term in (1.8) is modified to OL^{{\og{t— 

\xo\/v))/it-\xo\/v))-2). 

The proof of Theorem 2 is not discussed in the main body of this paper, since all 
of the needed information is contained in §4 and Appendix B of [3]. The main point 
is that Theorem 1 in fact establishes that for times \xq\/v + 1 <t < elogw, we have 

uix,t) = e-'*"'/2e**^/2e"'"NLSo(t - t2)[t(t;)sech(x)](x - xq - tv) 

+ e-^*^'/2g^i2/2g-^x•^,^LSo(t - ts) [r(f )sech(x)] (x + Xo + tv) 

where NLSo(t)0 denotes the free nonlinear flow according to (1.1). This is the start- 
ing point of the arguments provided in §4 and Appendix B of [3], which carry out 
an asymptotic (in time) description of the free nonlinear evolution of asechx, for a 
constant < a < 1. 

Although the main point of the present paper is to handle the difficulties involved 
in the case g < stemming from the presence of a linear eigenstate, some of the 
refinements we introduce (specifically, cubic correction terms in the interaction phase 
analysis) improve the result of [3] in the case g > 0. In fact, these refinements are 
simpler when carried out for g > directly, and we therefore write them out separately 
in that setting in §7. We summarize the results as: 

Theorem 3. In the case q > 0, the assumption (1.6) in Theorem 1 can he replaced by 
the less restrictive v > C(log |g|)^/^ + Ce_„(g)~, and the conclusion (1.7) holds with the 
first error term dropped (that is, C^L^dgl^w"^^^"^'^-') is dropped and only 0{v~^~^'^'^) is 
kept). Also, the conclusion of Theorem 2 holds with OL2(\q\3v~e(^~'^'^)) dropped from 
(1.9). 

Thus in the g > case we improve the L^ error from (9(t'~^/^+) to 0{v~^'^). It 
may be possible to to improve this error further to 0{v~'^~^) using an iterated integral 
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expansion of the error in the spirit of Sections 5 and 7, although a more detailed 
analysis than the one given there would be needed. 

We now outline the proof of Theorem 1, the main result of the paper, highlighting 
the modifications of the argument in [3] needed to address the case of g < 0. We will 
use the following terminology: the free linear evolution is according to the equation 
idtu + ^d^u = 0, the perturbed linear evolution is according to the equation idfU + 
\d1u — qSoU = 0, the free nonlinear evolution is according to the equation idtU + 
^d^u + |npn = 0, and the perturbed nonlinear evolution is according to the equation 
idtU + ^d^u — qSou + |npn = 0. 

The analysis breaks into three separate time intervals: Phase 1 (pre-interaction), 
Phase 2 (interaction), and Phase 3 (post-interaction). The analysis of Phase 2, dis- 
cussed in part earlier, is initially based on the principle that at high velocities, the 
time length of interaction is short ~ v~^~^, and thus the perturbed nonlinear flow is 
well-approximated by the perturbed linear flow. In [3], this was proved to hold for 
g > with a bound on the L^ discrepancy of size ~ w~2 + . In the case g < 0, we 
suffer some loss in the strength of the estimates due to the flow along the eigenstate 
|g|2e~l'sll^l^ and by directly following the approach of [3] the best error bound we 
could obtain is ~ |g|3t)~3+ -|- v~2 + . In the important regime |g| ~ v, this gives an 
error bound of size w~3 + , which does not suffice to carry through the Phase 3 post- 
interaction analysis discussed below. For this reason, we are forced to introduce a 
cubic correction term to the linear approximation analysis in Phase 2. The Strichartz 
based argument then shows that the L^ size of the difference between the solution 
and the linear flow plus cubic correction is of size ~ |g|3t;~6+ + t>~^+. However, since 
the cubic correction term is fairly explicit, we can do a direct analysis of it (not using 
the Strichartz estimates) and show that it is also of size v~^~^. Thus, in the end, we 
learn that the solution itself is approximated by the perturbed linear flow with error 

I |1 -1-1- -14- 

|g| 3v 6~^ -\- V . 

We then carry out the analysis of the perturbed linear evolution, as disscussed 
earlier, and show that by the end of the interaction phase, the solution is decomposed 
into a transmitted component (modulo a phase factor) 

(1.10) t{v)e"\edi{x - xo - t2v) 
and a reflected component (again modulo a phase factor) 

(1.11) r{v)e-'''''sech{x + Xo + t2v). 

In the post-interaction analysis, we aim to argue that the solution is well-approximated 
by the free nonlinear flow of (1.10) (that we denote Utr) plus the free nonlinear flow 
of (1.11) (that we denote Uj-d). It is at this stage that the most serious difficulties 
beyond those in [3] are encountered. The approach employed in [3] was to model the 
solution u as u = Utr + Wrcf + w, write the equation for w induced by the equations 
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for u, Utr, and u^cf, and bound lltf lU,, , ,^2 over unit-sized time intervals \ta,tb\ in 
terms of the initial size ||t(;(ta)||2,2 for that time interval. This was accomplished by 
using the Strichartz estimates. The Strichartz estimates provide a bound on a whole 
family of space-time norms ||iy||i'3 ^r where {q,r) are exponents satisfying an ad- 

missibilty condition - + j. = \- This family includes the norm L'^ t]-^x'i ^^^ other 
norms (such as Ln f.-iL^) are needed since they necessarily arise on the right-hand 
size of the estimates. From these estimates, we are able to conclude that the error at 
most doubles over unit-sized time intervals, and thus after ~ elogv time intervals, 
we have incurred at most an error of size v^. 

This strategy presents a problem for the case g < 0, since the linear eigenstate 
|g|i/2g-|g||x| is well-controlled in L?^ ^ ,L^ (of size ~ 1) but poorly controlled in L9^ ^ ,L^ 

(of size ~ \<l\^)- We thus opt to model the post-interaction solution as m = Wtr + 
Mrof + Mbd + w, where Wbd is the perturbed nonlinear evolution of the L^ projection of 
{u{ta} — Ure{{ta) — Mtr(ia)) outo the linear bound state |g|2e~l'^ll^'l. Then we can use 
nonlinear estimates based on mass conservation and energy conservation to control 
the growth of Ubd over the interval [ta,^b]- Then w{ta) is orthogonal to the linear 
eigenstate, and we can use the Strichartz estimates to control it over the interval 
[ta,^;,]. In the estimates, we take care to only evaluate Wbd in one of the norms 
controlled by mass or energy conservation. This argument is carried out in detail in 
§6. 

Acknowledgments. We would like to thank Maciej Zworski for helpful discussions 
during the preparation of this paper. The first author was supported in part by 
NSF grant DMS-0654436 and the second author was supported in part by an NSF 
postdoctoral fellowship. 



2. Scattering by a delta function 

Here we present some basic facts about scattering by a 5-function potential on the 
real line. Let g < and put 

The operator Hg is self-adjoint on the following domain: 

V{H,) = {fe H'iR \ {0}) : /'(0+) - /'(O") = 2g/(0)}, 

where /(O) means lima;^o/(3;) and /'(O^) means lim^^^oi f'i^)- This can be seen by 
verifying that the operators Hq±i are both symmetric and surjective on T>{Hq). We 
define special solutions, e±(x. A), to {Hg — X^/2)e± = 0, as follows 

(2.1) e±(x. A) = tqi\)e^'^'^xl + {e^'^^ + r,(A)e^*^^>'^ , 
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where tq and r^ are the transmission and reflection coefficients: 
(22) MA) = ^. r,(X) = -^^. 

They satisfy two equations, one standard (unitarity) and one due to the special struc- 
ture of the potential: 

(2.3) |t,(A)r + K(A)p = l, t,(A) = l + r,(A). 

Let P denote the L^-projection onto the eigenstate e*^'^'. Specifically, 

(2.4) P0(x) = |g|^/'e^l"l f \q\^/\'^^^(j){y) dy 

We have e~**^9P</)(x) = ea**"? P(j){x). Note that P(j) is defined for G L^., 1 < r < oo, 
and by the Holder inequality, 

(2.5) ||P0||l'-2 <c|g|^~^||0|U.i, l<ri,r2<oo. 

We use the representation of the propagator in terms of the generalized eigenfunc- 
tions - see the notes [8] covering scattering by compactly supported potentials. The 
resolvent 

R,{\) "^ {Hq - AV2)-^ , 
is given by 

RqWix, y) = -^-fj^ (e+(a^> A)e_(y, A)(a: - y)^ + e+{y, A)e_(a;, A)(x - y)^_) 

Using Stone's thoerem, this gives an explicit formula for the spectral projection, and 
hence the Schwartz kernel of the propagator: 

1 /"OO 

(2.6) e-'*^^ = -y e-^*^^/2(e+(x,A)M^ + e_(x,A)M^) rfA + e^^*'^^ 
We introduce the following notation for the dispersive part of e"**^'': 

The propagator for Hq is then described in the following 

Lemma 2.1. Suppose that cj) & L} and that supp C (— oo,0]. Then 

(2.7) e~''"''<P{x) = e'2^'1' P(f){x) + e-**^«((/. * Tq){x) xl 

+ (e-'*^"0(x) + e-**^°(0 * pq)i-x)) x^ 
where 



(2.8) Tq{x) = Soix) + pq{x), pq{x) = QC 



""x^ 
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Observe that we have, using a deformation of contour, 

p (A) =q [ e"(^-*^)cix = ^— f e'dx = rJX) 
Jo Q-^^ Jo 

f,(A) = l + r,(A)=t,(A). 

Observe also that HgR = RHg, where R(j){x) = 0(— x), so that the restriction on the 
support of is not a serious one, and the formula will allow us to estimate operator 
norms of Uq using e~**^*^ . Indeed, from the Hausdorff- Young inequality for e"**^" we 
conclude 



LP' 5 



(2.9) \\UMlv < 

where p G [2, cxd] and p' = -^. 
Proof. It is enough to show 

Uq{t)(t){x) = [e-'*^V(a;) + e-'*^°(0 * p,)(-x)]x'i + e-**^«(0 * r,)(x)a;^. 
From the definition of the propagator we have 

U,{t)^{x) = ^j y"e-**^'/2 (^e^(a;,A)^7(^ + e_(x,A)M^) <P{y)dydX, 

and so we must verify 

i- y" e-**^'/2 (e+{x, X)J^^^ij;X}<l>{y)dy + e_(x, \) j ^J^;^(l>{y)dy\ dX 

= [e-'*^V(a;) + e-**^°(0 * p,)(-x)]x^ + e-**^°(0 * rq)(x)x^. 
We compute first 

^^(^;^(t){y)dy= I e~'^y(^{y)dX + VJX) I e'^y<j){y)dX 

= 0(A) + r,(-A)0(-A), 

M^(/.(2/)% = t^lAj /" e^^V(i/)dA = t,(-A)<^(-A). 

J —oo 

We first verify the equation for positive x: 

e-^*^'/2/^e+(x,A) [^^{^)<P{y)dy + e.{x,X) f Jl^<t>{y)dy\ dX 



1^ 



1 r°'^ 





oo 



27r. 

- / e-^*^'/2 (^^^^^^ ^)^^^^) ^ ^^(_A)<^(_A)) + e_(x, A)t,(-A)0(-A) j dA 



SttJo 



1 



oo 






+ (e-*"" + r,(A)e*^^>,(-A)0(-A))(iA 



FAST SOLITON SCATTERING BY ATTRACTIVE DELTA IMPURITIES 



11 



At this stage we use tq{X)rg{-X) + rg{X)tg{-X) = 2Re (tq{X)rg{Xy\ = 0: 



1 

2'n 



— / e-**^'/2^,(A)e*^"0(A)rfA 



-itHo I 



r5*0)(x). 



The proof for negative x is similar, except that it uses rq{X)rg{—X) + tg{X)tg{—X) = 
|r.(A)P + |t,(A)p = l. D 

We have two simple applications of Lemma 2.1: the Strichartz estimate (Propo- 
sition 2.2) and the asymptotics of the linear flow e**^'' as f ^ oo (Proposition 2.3). 
We start with the Strichartz estimate, which will be used several times in the various 
approximation arguments of §3. 

Proposition 2.2. Suppose 

(2.10) idtu{x, t) + \dlu{x, t) - q5o{x)u{x, t) = f{x, t) , u{x, 0) = (f){x) . 



Let the indices p, r, p, r satisfy 



2 1 1 

2 < p,r < CO , I < p,f < 2 , — I — = - 

p r 2 



2 1 

- + - 
p r 



5 
2 



and fix a time T > 0. Then 



(2.11) 



I«IIl|'„_,jLs<c(||0|U2+t^||P0IU. 



4,,,.j + ^^ll^/IU;o,.,^s 



The constant c is independent of q and T. Moreover we can take f{x,t) = g{t)6o{x) 

4 



and, on the right-hand side, replace \\f\\j^p j^^ with\\g\\ 4 and replace Tp\\Pf\\ii ^r 



with Tp \q\^ r \\g\\^i 



"[0,T] 



[O.T] 



Proof. This will follow from 



(2.12) 



Ug{t)(l){x)+ I Ug{t-S)f{X,S)dS 





<c(||0||l2 



^10,T]^ 



TV Tf, 



To prove this estimate, we observe that the case g = is the standard Euclidean 
Strichartz estimate (see [7] and [3, Proposition 2.2]). The case q ^ ^ reduces to this 



case as follows. We write 



x_(j) and </)"'" (x) = Rx,(j) where again R(j){x) 



-X) 



Note that = 0" + /?(/>+, that Ug{t)R = RUg{t), and that {Rf) * g = R{f * Rg). 
Now, from Lemma 2.1, we have 

t/,0 =[f/o(t)0~ + f/o(t)(0~ * Pg)y^ + f/o(t)(0~ * rg)xl 

+ [f/o(t)0+ + f/o(t)(0+ * Rp,)\xl + f/o(t)(0+ * i?r,)a;°_. 
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We must now show that ||/^ * CTgllj;^!? j^f < c\\f\\jj, j^f, where a^ is either r^ or pg. 



"[O.T]-' 



This follows from applying Young's inequality to the spatial integral. 

This completes the proof of (2.12). To obtain (2.11), we observe that ||P(/)(x)||x,p ^r 



ri/p||P0(x)|U. and j'Pf{x,s)ds 



TV rr 



<TVp||P/(x,s) 



"[O.T]-' 



[0,T]^ 

• . The first is im- 



mediate, and the second follows from the generalized Minkowski inequality: 



P/(x, s)ds 



< 






\\Pf{x,s)\\Lr^ds 



[0,T] 



<T'/P I \\Pf{x,s)\\Lrds. 



D 



We now turn to the large velocity asymptotics of the linear flow e" 



-itH„ 



Proposition 2.3. Let 9{x) he a smooth function hounded, together with all of its 
derivatives, on R. Let cf) G 5(M), v > 0, and suppose supp[^(a;)0(x — Xq)] C (— oo, 0]. 
Then for 2|xo|/f <t<l, 
(2.13) 

-itHo l^—ixv 



e-^*^'[e^^Xx - xo)] = t(t;)e-**^°[e"Xx - Xq)] + r{v)e- —'\e 



-itHo To*^*'^ 

+ e(x, t) 



-X - Xo)] 



where, for any k > 0, 



|e(x,t)|U. < kdMxmx-xo)]\\L^^ + j^\\{xmx)\\ 

U \LUj 



m 



4||(1 - e{x))<P{x - x^)\\li + ||P[e"^^(x)0(x - Xq) 



\Ll 



Remark 2.1. In §3, Proposition 2.3 will be applied with 9{x) a smooth cutoff to 
X < 0, and 0(x) = sech(x) with xq = —v'^ <C 0. 



Before proving Proposition 2.3, we need the following 
Lemma 2.4. Let tp G iS(M) with supp ^ C (— cx3,0]. Then 



(2.14) t/,(t)[e^^'^^(x)](x) 

where 

uniformly in t. 

Proof. We apply (2.7) with > 



-itHo [0*2:1) 



ip{x)]{x)x^_ + t{y)e 



-itHo [0*2:1) 



ip{x)]{x)x^. 



+ r(i;)e-**^ne-*^^ 



il){—x)]{x)x^_ + e(x, t). 



|e(x,t)||L2 < -\\d^i)\\L2 

V 



[X) 



e^^'"il){x) to flnd that 



e(x, t) 



-itHo 



[(p * {pq - rq{v)6o){-x)]x^_ + € 



-itHo [ 



*(^<? 



tg(w)(5o)(x)]x° 



We pass to the Fourier transform using Plancherel's theorem: 



\e-''"°ip\\L2 = c\\ip\\L2. 
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So that it suffices to verify 

||^(A-t;)(r,(A)-r,(t;))|U. <^||A^(A)|U., 

||^(A-^;)(t^(A)-t,(^;))|U. <^||A^(A)|U.. 
We write out now 

[iX — q)[iv — q) [iX — q)[iv — q) 

/ M , /, N / N , |A — Wl lA — f I 

\r,w - .,(«)! ^ i*,(A, - *,(.,,! ^ ^,,^;^,;„.^^, < ^- 

We plug tliis in to obtain our estimate: 

||^(A-tO(r,(A)-r,(i;))|U. = ||^(A-tO(t,(A) -t,(i;))|U. < ^||A^(A)|U. < ^WOMl-- 

U 
Now we turn to the proof of Proposition 2.3. 
Proof. We will prove (2.13) by showing that 

Uqit) [e^^Xx - xo)] = t(w)e-^*^o [e^^Xx - Xq)] + r(w)e-**^o [e-™^0(-x - Xq)] + e(x, t), 
where, for any A; > 0, 

||e>,t)|U. <-||9.[^(x)0(x-xo)]|U| + ^||(x)V(x)|k.+4||(l-e(x))0(x-xo)|Ui^ 

We use (2.14) with 'il){x) = 9{x)cf){x — Xq), which gives (using ei for the error term 
arising from the lemma), 

f/g(t)[e"-''^(x)0(x - Xo)](x) =e-'*^«[e""^(x)0(x - Xo)](x)x^ 

+ t(i;)e~**^°[e"^^(x)0(x - Xo)](x)x^ 
+ r(i;)e-**^°[e-"^^(x)0(-x - Xo)](x)x'i 
+ ei(x,t). 

We rewrite this equation with 9 omitted at the cost of an additional error term: 

f/g(t)[e"Xx - Xo)](x) =e-^*^«[e*^Xx - Xo)](x)x^ 

+ t(t;)e-'*^°[e'"Xx - xo)](x)x^ 

+ r(^;)e-**^«[e-"X-x - xo)](x)x'i + ei(x,t) + e2(x,t). 

Using the notation /(x) = e"'*'(l — ^(x))0(x — xq), this error term is given by 

e2(x,t) = Uq{t)f{x)-e-''^'>f{x)x'l-t{v)e-''^''f{x)xl-r{v)e~''"°f{-x)x^_. 
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Recall that e ^*^° is unitary, and that \\Uq{t)\\Li^Li = 1- This gives us a bound on 
62: 

\\e2ix,t)\m < 4||/IUi = 4||(1 - 9ix))4>ix - xo)^. 
We have a bound on ei from the lemma, so combining everything we have, we see 
that to prove the proposition it remains only to prove 

||e-^*^°[e^^Xx - xo)]{x)x'_\\l2 + ||t(tOe-^*^°[e^^^^(x)</)(x - xo)](x)x°_|U| 
+ ||r(tOe-*^ne""->(-x-Xo)](x)x°JU. < ^|Kx)V(x)||^., 

for every A; > 0. However, because \t{v)\ < 1 and 1^(17)1 < 1, and because e~^^^°R = 
Re~^^^° , where R(j){x) = 0(— x), it suffices to prove 

(2.15) ||e-*^°r>(x-Xo)](a;)x°_|U. < ^||(a;)V(x)|k.. 

We expand as follows, first using the definition of the propagator: 



-itHo \ ixv 



2n 



X - xo)](x) = — / e"^-**^'/2 / e-'^ye'y-(j)(y _ XQ)dyd\. 



Here we change variables y ^^ y + Xq io simplify the Fourier transform. 



27r 



,ixA-itAV2gixo(^-A)^(^^ _ ^^)^^_ 



Here we change variables \^^ \ + v. 
1 



27r 



ixv I ix\—it(\+v)/2—ixo\ 



4>{\)d\ 



= _^ixv^-itvy2 / e*^(^-^o-*'')e-^*^'/20(A)rfA. 
27r J 

After k integrations by parts in A, this becomes 

1 pixv—itv'^j2 



x(x-xo-tv)Qk L-^t^y^^(x)'\ dx. 



27r {i{x — Xq — tv)Y 

By assumption we have 2\xq\/v < t < 1, so that — xq — tv < Q and |a; — xo — tv\ > 
I — xq — tv\ > tv/2 when x < 0. Hence 

Cfc 



(x-Xo)](x)x_ 2 < 



Ck 



iX(x-xo-tv)Qk fe~''^'/^^{X)'] dX 



Ll 



{tvf 

k 



X 



r2 



\m- 



J=0 



Since t < 1, (2.15) follows. 



D 
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3. SOLITON SCATTERING 

In this section, we outline the proof of Theorem 1, the details of which are executed 
in §4-6. We recall the notation for operators from §2 and introduce short hand 
notation for the nonlinear flows: 

• Hq = —\d1. The flow e~**''^° is termed the "free linear flow". 

• Hq = —\d1 + g5o(x). The flow e"**^"? is termed the "perturbed linear flow". 

We also use Uq{t) = e"**^' — ea**^ P, the propagator corresponding to the 
continuous part of the spectrum of Hg. 

• NLSq(t)0, termed the "perturbed nonlinear flow" is the evolution of initial 
data 0(x) according to the equation idfU + ^d^u — qSo{x)u + \u\'^u = 0. 

• NLSo(t)0, termed the "free nonlinear flow" is the evolution of initial data 0(x) 
according to the equation idth + ^d^h + \h\'^h = 0. 

From §1 and the statement of Theorem 1, we recall the form of the initial condition, 
Uo{x) = e"''sech(x — Xq), w ^ 1, Xq < — f^ where < e < 1 is fixed, and put 
u{x,t) = NLSg{t)uo{x). We begin by outlining the scheme, and will then supply the 
details. In this section, the O notation always means L^ difference, uniformly on the 
time interval specified, and up to a multiplicative factor that is independent of q, v, 
and e. 

Phase 1 (Pre-interaction). Consider < t < ti, where ti = \xq\/v — t)"*^^""^) so 
that Xq + vti = —v^. The soliton has not yet encountered the delta obstacle and 
propagates according to the free nonlinear flow. Indeed, there exists a small absolute 
constant c such that {q)^e~'^'"^ < c implies 

(3.1) u{x,t) = e-'*"'/2git/2g*x.gg^j^^^ _^^_ ^^^ ^ ^^|^|3^^^i^_„e^^ < t < ti. 

This is deduced as a consequence of Lemma 4.1 in §4 below. 

Phase 2 (Interaction). Let ^2 = ti + v^^"^, so that Xq + t2V = v^ , and consider 
ti < t < t2. The incident soliton, beginning at position — f^, encounters the delta 
obstacle and splits into a transmitted component and a reflected component, which 
by t = ^2, are concentrated at positions f ^ and — f '^, respectively. More precisely, there 
exists a small absolute constant c such that if i;(q')^ ( e 2 + e~'^'"'' 1 +w~3(^~^)|g| 3 < c, 
then at the conclusion of this phase (at t = ^2), 

(3.2) m(x, h) = t(^;)e-**2^'/2git2/2gix,;gg^j^^^ _ ^^ _ ^^^^ 

+ r(^;)e-**2^'/2git2/2g-ix«gg^^(^ + Xq + vt2) 

+ C(tri(i-^)|g|^) + C(tr(i-^)). 
This is proved as a consequence of Lemmas 5.1, 5.2, and 5.3 in §5. 
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Phase 3 (Post-interaction). Let ^3 = ^2 + £logf, and consider [t2;^3]- Suppose 
|g|i0|/g,^3^,-i4(i-2e) _|_ ^|/g^3 /g ^ + e"'^^'^ j < c and {q)v~"' < c^^n-, where c is a small 
absolute constant and c^^n is a small constant, dependent only on e and n, which goes 
to zero as £ ^ or n ^ 00. The transmitted and reflected waves essentially do not 
encounter the delta potential and propagate according the free nonlinear flow, 

(3.3) n(x, t) = e-'*^''/2e**2/2e«^NLSo(t - ta) [t{v)sech{x)] {x - Xq - tv) 

+ e-'*"'/2e'*^/2e-'"''NLSo(t - t2)[r(f)sech(x)](x + xo + tv) 

+ Oiv~^'-'^) + 0(|g|^^;-i(i-2^)), h<t< ts. 

Now we turn to the details. 

4. Phase 1 

Let Ui{x,t) = NLSo(t)uo(x) and note that 

Mi(x,t) = e-'''''^^e''/^e'''''sech{x - xo - tv) 

Recall that ti = |a:o|/t' — v~^^~'^\ so that at the conclusion of Phase 1 (when t = ti), 
the position of the soliton is xq + vti = —v'^. Recall that u{x,t) = 'NLSq{t)uo{x) and 
let w = u — ui. We will need the following perturbation lemma. 

Lemma 4.1. Ifta < tb < ti, tb — ta < ci < I, and 

{q)^\\w{x,ta)\\Ll + \q\Hq)^\\ui{0,t)\\Li ^ , < C2|g|"^ 

then 

\\w\\li Lr<C3({q)^\w{x,ta)\\Ll + \q\Hq)^\ui{0,t)\\Li 

where {p, r) is either (00, 2) or (4, cxd), and the constants ci, C2, and C3 are independent 
of the parameters e, f , and q. 

Before proving this lemma, we show how the phase 1 estimate follows from it. Let 
/c > be the integer such that kci < ti < [k + l)ci. (Note that k = ii the soliton 
starts within a distance v of the origin, i.e. —v — v^ < Xq < —v"^, and the inductive 
analysis below is skipped.) Apply Lemma 4.1 with ta = 0, tb = ci to obtain (since 
u;(-,0) = 0) 

IklU-^^jLi < C3|g|^(g)^||ni(0,t)||L^_^^j < C3|g|i(g)^sech(xo + wci). 

Apply Lemma 4.1 again with ta = ci, tb = 1c\ to obtain 

lklU-_^^^,L2 < c3((g)^||w(-,ci)|U| + |g|^(g)^lhi(o,t)|U-^^^j) 

< Cglgl^ (g)sech(xo + vc\) + Cglgj^ (g)2sech(xo + Ivcx). 
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We continue inductively up to step k, and then collect all k estimates to obtain the 
following bound on the time interval [0, kci]: 

k 

IklU-^^^jLi <C3|g|5^c3"^(g)^sech(a;o+jwci). 
i=i 

We use here the estimate secha < 2e~'"': 

fc-i 



<C3^|g|^(g)^e-«+^-5^C3-^'(g)- 



'tgJI'Cl^ 



j=0 



We introduce here the assumption that Cg ^{q) 23"^^ > 2, which allows us to estimate 
the geometric series by twice its final term, giving 

Finally, applying Lemma 4.1 on [A;ci,ti], 

|g| 2(g) 2 6^0+^1^'= + \q\2{q)2sech{xo + tiv)j < c\q\^{q)2e-''\ 

where we have used Xq + Civk < Xq + tiv = —v^. The constant c here is still 
independent of g, v, and e. As a consequence, (3.1) follows. So long as this last 
line is bounded by c{q)^ the repeated applications of Lemma 4.1 are justified. 
Now we prove Lemma 4.1: 

Proof. We begin by writing a differential equation for w in terms of ui. 

idtw + -dlw — qSo{x)w = ui\ui\'^ — {w + ui)\w + wip + qSo{x)ui 

= —w\w\ — 2ui\w\ —Uiw — 2\ui\ w — u^w + qoo{x)Ui. 

We rewrite this as an integral equation in terms of the perturbed linear propagator 
^-iHqt _ u^[i^ _)_ g«*9 /2p^ regarding the right hand side as a forcing term. 

w{x, t) = {U,{t - to) + e^(*-*'')^^/2p) nj{x, Q 

+ / Ug{t — s) [—wlwl"^ — 2ui\w\'^ — Uiw"^ — 2\ui\'^w — ufw + q5o{x)ui^ ds 

Jta 

Jta 

= I + II + III. 

We define \\w\\x = ||'W^||l°° l^ + ||w|U4 ^oo, and then proceed to estimate the X 
norm of the right hand side term by term. In what follows, (p, r) denotes either (00, 2) 
or (4, 00). 
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I. We observe from the Strichartz estimate that 

\\Uq{t - ta)w\\x < c\\w{x,ta)\\Ll- 

Next, using the Holder estimate for P, we have 

Taken together, and using tf, — ta < ci, these bounds can be written as 

II. For the terms involving Ug we will use the Strichartz estimate, which tells us that 



/ Uq{t- s)f{x,s)ds 

Jta 



X 



-'''-"'^f;..,^s' 



whenever 2/p+ 1/r = 5/2. The first term, w;|tf;p, is cubic, and we will use p = 1 and 
f = 2: 



Uq{t — s)w\w\'^{x, s)ds 



Jta X 

We first pass to the L^ norm for two of the factors. 



< Cllw II rl T2. 



< c\\\H\Ll{t)\\w^L^{t)\\ , 



[ta.tft] 



And then pass to the L^ ^ , norm for the other factor. 

< c||w||l^ L2||w ||rl Tx = c\\w\\l°° L'^\\w\\t2 roo- 

Finally we use the boundedness of th — ta to pass to the X norm. 

<c{tb-ta)^\\w\W. 

The quadratic and linear terms follow the same pattern. Observe that the distinction 
between w and wJ, and between ui and Mi, does not play a role here: 

' Uqit- S)ui\w\^{x,s)ds <c{tb-ta)'^\\ui\\x\\w\\\<c{ti,-ta)'^\\w\\\ 



Uq{t — s)w\ui\^{x, s)ds 



X 



< c{tb-ta)2\\w\\L6 



,LS- 



For the delta term we have no flexibility in the choice of exponents. 

' Uq{t — s)5o{x)ui{0, s)ds < c\q\\\ui{0,t)\\^4/3 . 



"['a.tj,] 



III. The P terms we will similarly estimate one by one: 
/ gi(i-s)9V2p(^ _ ^|^|2 _2ui\w\^ -Ww^ - 2\ui\^w - ulw + qSoix)ui)ds 



ft™. til a; 
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< c(\\P\w\^\\rl rr + \\PUl\w\'^\\rl jr 

+ \\P\m\M\Ll ^ rj + |g|||P5o(x)ni(0,t)|Ui .. 

Here we used a generalized Minkowski inequality to pass the norm through the inte- 
gral, just as we did in the proof of the Strichartz estimate. Note that the constant 
c in the second line depends on p, but since p only takes two different values we will 
not make this dependence explicit in our notation. For the delta term as before we 
have no flexibility in the choice of U^ norm: 

|g|||P5o(xK(0,t)|Li ,,Ls<c|gP~'/lki(0,t)|Li <c|g|t(g)^||«i(0,t)||M ,- 

For the cubic term in w we proceed using the same Holder estimate for P. Here we 
use ri = 2, giving this term a factor no worse than (g)2: 

ll^kl'lki ^,Li<c{q)Hw^L] ,M = c{q)Hw\\\. 

The last step is the same as that in the w^ term in II above. We now estimate the 
quadratic and linear terms in w. We have 

llPMlW^IIri jT < C\\w'^\\ji rr = cllwll f2 r2r < Cllwll^ 

\\Pu^w\\t1 tt < c\\w\\t1 jr < cit}) —ta)^\\w\\x- 

When r = oo this last step is achieved by passing from Ll ^ , to Ln ^ , and using 
the boundedness of the time interval. When r = 2 we interpolate using Holder's 
inequality between L^ and L^. 

Having estimated each of the terms individually, we combine our results. We use 
tb — ta < Ci < 1 io pass from lower norms in time to higher ones, only tracking the 
power of (tfe — to) for the linear term in w. 

\\w\\x < c({q)^\w{x,ta)\\Ll + |g|t(g)^||ni(0,t)|Ui^^_^^^ 



+ {tb-taf'^\\w\\x + \\w\\\+\q\ 



1/3 I 



w\ 



X 



We now take ci sufficiently small (recall that t;, — t^ < ci) so that the linear term in 
||w||x can be absorbed into the left hand side: 

\\w\\x < c [{q)--\\w{x,ta)\\Li + |g|^(g)^||ni(0,t)|Ui^^_^^^ + ||w;||^ + \q\^'^\\w\\],^ , 

with a slightly worse leading constant c. We rewrite this inequality schematically 
using X = \\w\\x'- 

<A-x + Bx^ + Cx^. 

We now consider ||w|U(t') for ||w|U(i') = IklU- ^,jL2 + IkllL*^^,^^- for t' ^ [ta,h]- 
This is a continuous function oft', and, for each t' G [ta,^;,], ||w||x(i') obeys the above 
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inequality. Therefore if we find a positive value Xq for which the inequality does not 
hold, we will be able to conclude that ||w|U{t') < xq for every f G [ta,tb], and hence 
also that \\w\\x < xq. 

We will use Xq = 2A, and arrange A, B and C so that this gives a negative right 
hand side. In fact, we have 

A-2A + ABA^ + SCA^ = -A + A{4BA + 8CA^). 

To make this negative we impose 45^4 < j and SCA^ < i. We thus obtain x < 2 A, 
or, in the language of ||u^||x5 

\\w\\x <c(^{q)^\\w{x,ta)\\^ + \q\^q)^\\ul{0,t)\\^^^^^^^ , 
provided that {q)^\w{x,ta)\\Li + |g|2 (g)^ ||ni(0, t)||ii < c|g|~^/6. D 



5. Phase 2 

We begin with a succession of three lemmas stating that the free nonlinear flow 
is approximated by the free linear flow, and that the perturbed nonlinear flow is 
approximated by the perturbed linear flow. The first lemma states that the nonlinear 
flows are well approximated by the corresponding linear flows, the second gives a 
better approximation by adding a cubic correction term, and the third shows that the 
improvement is retained even if the cubic term is omitted. In other words, we 'add and 
subtract' the cubic correction. Our estimates are consequences of the corresponding 
Strichartz estimates (Proposition 2.2). Crucially, the hypotheses and estimates of 
this lemma depend only on the L^ norm of the initial data 0. Below, the lemmas are 
applied with 0(x) = u{x,ti), and ||m(x, ti)||L2 = ||mo||l2 is independent off; thus v 
does not enter adversely into the analysis. We first state the lemmas and show how 
they are applied, deferring the proofs to the end of the section. 

Lemma 5.1. Let cf) e L"^ andO<h. If {t^^ + tl^^\q\^^^) < Ci(||0|U2 +ty^||P0||i6)-2, 
then 

(5.1) ||NLS,(t)0 - e-''^''<Ph.^^^^,n < c,t]!\\ + t;/^|gp/^)(||0|U2 + tT\WU)\ 

•where P and R satisfy p + -^ = |, and Ci and C2 depend only on the constant appearing 
in the Strichartz estimates. We alert the reader that in our notation P is used both 
as a Strichartz exponent and as the hound state projection (2.4). 
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Lemma 5.2. Under the same hypotheses as the previous lemma, 
(5.2) 



|NLSo(t)0 — oil roc r2 < C 



1/6: 1 1/3 



L2+t 



1/6 
h 



L6 



tl(}^tT\^ 



+ tf fl + tf|g|V3) (11011,.+ t;/«||P0 



IlS 



1^2 



+i.(i+ifkr/')(ii0iiL2+tfiiP0ii 



L«fii0iiL 



where 



g(t) = e"'*^«(/)+ / e 
'o 



-i{t-s)Hq\ -isHqi\2 -isHq 



Lemma 5.3. For ti < t2 and = u{x,ti), we have 
(5.3) 



t2 



-i{t2-s)Hq\ -isHql\2 -isHq 



LI 



<c 



(t2 - ti) + (t2 - ti)^ ((g)^|g|^e-^^ + (g)>|V2^^ + (g)^|g|^e-3^^) 
where c is independent of the parameters of the problem. 

In order to apply these lemmas, we need to estimate ||P(;A||j,6. As before, let 0i(x) 



-itv /2 pit/2 pixv 



e**/ e*^''sech(x — xq — vti) and 02 



L. It suffices to estimate ||-P0i||l6 and 



II02IU? <c(g)2|g|2e ''' 



||P02|U6. By (3.1), 

and thus, by (2.5) 

\\P<p2U<c{q)^q\le-''\ 
On the other hand, a direct computation gives 



|-P0i(a^)| < Igle'''^' / e'^^^^sech.{y — Xq — vti)dy 



\q\e' 



q\x\ 



CQ+TJt]^ 



e'"^'sech(y — xq — vti)dy 
+ / e^'^'sech(|/ — xo — fti)(iy 



In the ffist integral we use the fact that e'^'^l is uniformly small in the region of 
integration. For the second we use the fact that e'^'^' is bounded by 1 and that 
secha < 2e~". 



|P0i(x)| < c|g|e^l^l /'e¥(^o+^'ti) + ^\[xo+vH)\ < ^^^^^q\x 



|g|i,= 



c|g|e^' ' e 2 + e 2 
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This implies that 

Combining, 

(5.4) 



||^0i||l6 < C|g|6 6 2 +6 2 



hhl 



r0l|LB<c(g)i|g|t(6-^^+6t-^). 

Set t2 = ti + 2v^~^ and 0(x) = m(x, ti). We now give an interpretation of our three 
lemmas under the assumption v{q)^ (e 2 + 6^~ ) + f~3(^~^)|g|3 < 1. This makes 
(5.2) into 

\\U{XM) - g{t2)\\Ll < c(v-' + V^l'^^-'W^y 

and (5.3) into 



9{t2 



<t2~t^)H,u{x,ti) 



< cv-^^-'\ 



Ll 



Overall this amounts to 

u(-,t2) = NLSg(t2-ti)[n(-,ti)] 

= 6-^(*^-*^)^'[u(-,ti)]+0(i;-i(^-^)|g|^) + 0(tr(i-^)). 

By combining this with (3.1) we find that under our new assumption the errors from 
this phase are strictly larger, giving 



(5.5) 



M(-,t2) 



e-i^^-V2^ihl2^-i{t,-t,)H,^^ixv^^^^(^^ _^^_ ^^^^ 



+ 0(i;-i(i-^)|g|i) + C(t;-(i-^)). 



By Proposition 2.3 with 9{x) = 1 for a; < — 1 and 6{x) = for x > 0, 
sech(x), and Xq replaced by xq + tiv, 

e-'(*^-*^)^'[6'"''sech(x -xo- vh)]{x) 

= t(i;)6-^(*^-*^)^°[6"^sech(x - Xo - vti)]{x) 



(X 



(5.6) 



+ r(^;)6-^(*2"*i)^°[6-*^^sech(x + Xq + vti)]{x) 



+ 0{v-^] 



where we have used the assumption that e~^^ < v~^. Now we use (7.2) and (7.3) 
to approximate 6"**^" by NLSq, picking up an error of 0{v^~'^). By combining with 
(5.5) and (5.6), we obtain 

u{-,t) = t(t;)6-^*^'''/2giti/2NLSo(t2 - ti)[6^^''sech(x - Xq - vh)]{x) 

+ r(t;)6-'*^'''/2g*ti/2^LSo(t2 - ti)[6-"''sech(x + xq + vh)]{x) 
+ 0(i;-i(i-^)|g|^) + C(tr(i-^)). 
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By noting that 

NLSo(t2 - ti)[e*^"sech(x - xq - hv)] 

and 

NLSo(t2 - ti)[e-*^''sech(x + xq + hv)] 

we obtain (3.2). 

Now we prove Lemma 5.1: 

Proof. Let h{t) = NLSg(t)0, so that 

idth + -dlh - q5o{x)h + \h\'^h = 0. 

' 4 r oo 



Where in Phase 1 we used L^L"^ as an auxihary Strichartz norm, here we will use 

;e the nol 
Strichartz estimate 



L^L^.. We introduce the notation \\h\\x' = ||^||l°° lI + II^IIl^ ^6, and apply the 



once with (p, r) = (oo,2) and (p, f) = (6/5,6/5), and once with (p, r) = (6,6) 
and (p, f) = (6/5,6/5). We observe that Holder's inequality implies that ||/||ip < 
\Lvi\\f\\LV2\\f\\LV3 provided ^ = ^^+^ + ^^. This gives us 







\h\\. 


L8/5 ^18/5 


< \h\l. 


igi^i^, 


..^l ^ ^''^ 


'^\h\\,. 


We also have 






















\P\h\^h\ 


^[O.ttl^- 




< .1/2 

LI ^h 


Wx'. 








^1/6, 


\P\h\^h 


^[O.tfcl^^ 


< ctr , '' 










yielding 


















\h X' 


<c( 


L2 + 


\P4>W 


+tf |P0 


L^+t'/- 


'h\\, 


+ t 


I'W'" 



We then use the fact that P is a projection on L^ to write 

\\h\\x' < c{U\W + tf IIP0IU6 + tl''\\h\\\, + tl/\\''^\\h\\\, 
Using, as in Phase 1, the continuity of ||/i||x'(tfc)5 we conclude that 

\\h\\x> <2c{U\y+tl"\\P4>\\,.l 
so long as %c\tl'^ + tl'\\^l^){U\y + C\\P<P\yf < I. 
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We now apply the Strichartz estimate to u{t) = h(t) — e~**^90, observing that the 
initial condition is zero and the effective forcing term —\h\^h, to get 

(5.7) ||Mt)-e-^*^^0|L.^^^,«<c|||/.r/.||,a/s^^,a/5+t;/^|^^ 






Now we prove Lemma 5.2: 
Proof. A direct calculation shows 

h{t)-g{t) = [ e-'(*"^)^' {\h{s)\^h{s) - |e"*^^'0pe-*^^'0) ds. 
Jo 

The Strichartz estimate gives us in this case 



D 



isHa ±\2 „—isH, 



\h-g\\L^^^r.<\\\h\'h-\e-'^^^<P\'e 



6/5 f6/5 






\p{\h\^h-\e-'^^^<P\^e-'^''^cp)U}..,Li 



"[O.ti,]' 



I + II. 



We introduce the notation w{t) = h(t) — e ^'^^''(p, and use this to rewrite our difference 
of cubes: 

\h\^h - \e-''"''(j)\^e-''^''(f) = w\w\'^ + 2e-''^''(f)\w\^ + e''^''4>w'^ 

isHn 1\2„,. I I „—isHa l\^ rr-r 



+ 2 |e-*^^'0 !'«;+( 



w. 



We proceed term by term, using Holder estimates similar to the ones in the previous 
lemma and in Phase 1. Our goal is to obtain Strichartz norms of to, so that we can 
apply Lemma 5.1. 

I. We have, for the cubic term. 



-') 
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For the first inequality we used Holder, and for the second Lemma 5.1. Next we treat 
the quadratic and linear terms using the same strategy (observe that as before we 
ignore complex conjugates): 









< Ctf (1 + tf |g| V3)(||0||,, + ty6||p^||^^)6||^|| 



L2. 









^1/6 I 



<ct,(||0|U.+tri|P0|U6)^||0||i.. 

II. In this case we have 

I2„.,ll . / „+l/2|L„||3 

(.i/6|iD^ii ,.^9 









i2+t^/''||P0|U6) 



r6 r6 



\Pw\e 



-isHa Al2l 



1/2 I 



L2, 



rl r2 < Ct, kf r6 r6 6 

-^[o.t^]^^ - ^ II ll^lo.t^i-^-ll 



L2 



,1/6 



+tr\w\\L^r\mh- 



Putting all this together, we see that 



l/l — oil TOO 7-2 < C 






L2 



+ t,(l + t;/^|g|V^) (11011,2 +tf 



Le)l0|li2 



Finally we prove Lemma 5.3: 



D 



Proof. We write (j){x) = 0i(x) + 4>2{x), where 0i(a;) = e~**^*' /^e**^/^e*^''sech(x — Xq — 
uti), and estimate individually the eight resulting terms. We know that for large v, 
02 is exponentially small in L^ norm from Lemma 4.1. This makes the term which is 
cubic in 0i the largest, and we treat this one first. 



I. We claim 



/t? e-^(*2-^)H,|g-i.//,0^|2g-isH,0^^g < ^j^^^ _ ^^y 
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We begin with a direct computation 



(5.8) 


Jtl 


LI 


(5.9) 


< c(t2 - ti) |e-*^^^</)i 2e-^^^'0i 


C Vl TOO r2 



It remains to show that this last norm is bounded by a constant. We use (2.7) to ex- 
press e"**^' in terms of e"**^", recalhng the formula here for the reader's convenience: 



e 



i(x) = [e-^*^Vi(x) + e-^*^°(0i * p,){-x)y_ 



This formula is only valid for functions suported in the negative half-line, but this 
will not cause serious difficulty and we ignore the problem for now. We first evaluate 
this expression with e^^'"il){x) in place of 0i(x). Here il){x) = sech(x — xq — vt), and the 
other phase factors do not affect the norm. The first term uses the Galilean invariance 
of e~^^^° directly: 

For the second and third terms we use in addition the fact that e~^^^° is a convolution 
operator, and convolution is associative: 

e-^*^o(e^^^^*Pg)(-x)x^ = [(e-^*^Oe*^» * pg] (-x)x^ 

= e-^*^'/2 [(e*^''e-^*^°7/;(x - vt) * pg] {-x)x^_ 

[{e-''"yix - vt) * (e-*"»] (-x)x'i 



^-itv'^/2^-ixv If^-itHo^ 



e-^*^o(e"'> * Tg){x)xl = e-'^^'/^e"" [{e-''"°ij{x - vt)) * (e'^V^)] {x)xl 



The final term we leave as it is, so that we have 

(5.10) e-'^^''e''"'ij{x) =e-**'''/2gix« [e-^*^V(a: - vt) 

(5.11) + [{e-''^'>^{x - vt) * (e-*^»] i-x)]x°_ 

+ e-^^'^'/Sgix,; [(g-it^o^^^, _ ^^)) ^ (e""Vg)] {x)xl 

(5.12) +e5**'''Pe™^7^(x). 

Before proceeding to the estimate of (5.8) we introduce the following notation: 

/_(x) = /(x)x°_, /+(x) = /(-x)xO„ 

so that / = /?/+ + /_, where Rf{x) = /(— x), and (5.10) will be applicable to f±. 
We then write 



\e-''^^(j)i\^e-''"''(f)i = R{\e~''^^(pi+\'^e-''^^(j)i+) + ■■■ + \e-''^^(j)i_\'^e-''^^c 



>i-, 
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with a total of 8 terms. After applying (5.10) three times to each of them, we will 
have 8 ■ 4'^ terms. We will estimate these terms in groups. We observe first that the 
distinction between 0i+ and 0i_ will not play a role, and neither will the presence or 
absence of R. We accordingly write ^ for sech(a; — Xq — vti)± and omit R when it 
appears. In what follows aq denotes either e~^^'"pq, e~'^^'"Tq, or 6o, each of which has 
L^ norm 1. 



-isHo 



1p{x — VS)) * CTq 



-isHi 



°%l){x — vs)) * a^^ 






< c 



{e-'''"'ij{x - vs)) * aJ {e-''''y{x - vs)) * cr, 



isHo„ 



r oo r 2 



Now we pass from L^ to L°° for two of the factors, and use Young's inequality: ||/ * 
g\\p < ll/llplls'lli) once for p = 2 and twice for p = oo, and then use the Gagliardo- 



Nirenberg-Sobolev inequality which states that the L°° norm is controlled by the H^ 
norm. Because the H^ norm is preserved by e~**'^°, we are home free: 



< c (e- 



-isHo 



ip{x — vs)) * a, 



iWl 



f 2 

[tl,t2l^ 



\\{e-''"°i){x-vs))*(T, 



Q r oo r oo 

"^[tl,t2l^ 



< C||sech(x)||roo r2 ||sech(x)|| roo rrl < c. 



Terms where one or more of the factors of (e ^'^^°il)[x—vs)) are replaced by ea***? Pe^^ijji^) 
are treated in the same way. We have 



e¥t'i^ Pe''"'^{x) 



< ||^(a:)||roo rp, 

iff . ,iS ~ ^"^[tl,t2l^-' 

[«l.t2l 



where p is either 2 or cxd. 

II. For the other terms the phases will play no role, because we will use Strichartz 
estimates. The smallness will come more from the smallness of 02 than from the 
brevity of the time interval. 



ft2 

Jtl 



^ III ~isH„ 1 1 2 ~isH„ 

< c\\\e '01 e '( 



Li 



'2lli[ti,t2]-C'i- 



We have used the Strichartz estimate with (p, f) = (1, 2), and combined the resulting 
terms using (2.5). We use Holder's inequality so as to put ourselves in a position to 
reapply the Strichartz estimate. 



< c(t2 -ti)2||e 



-isHn 



'^"4..2l^S°ll^ 



-isHn 



'2ll^[?,.,,^i 



< C(t2 - tl)M II01IU2 + (t2 - tl)3||P0 



l\\L°^ 



2\\li + \\PUli). 
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Once again we use (2.5) to combine terms, this time with a penalty in |g|. 

<c(t2-ti)^(g)||0i||i.||02||Li 
< c(t2 -ii)2(g)2|g|2e"'''. 
Similarly we find that 



Ju 



2\„\2,„-2v'' 



<c{t2-h)--{qY\qYe 



LI 



i2 



^-i(t^.s)H^^^-isH,^^^2^-isH^^^^^ 



< c(t2 -ti)2{q)2\q\2e' 



-3v^ 



Ll 



u 



6. Phase 3 

Let ^3 = ^2 + elogf. Label 

Utr{x,t) = e-'*'''/2gif2/2gix.NLSo(t - t2)[t(i^)sech(x)](x - Xq - tv) 

for the transmitted (right-traveling) component and 

u,,i{x,t) = e~**'''/'e'*^/'e-*"''NLSo(t - t2)[r{v)sech{x)]{x + Xq + tv) 

for the refiected (left-traveling) component. By Appendix A from [3], for each fc G N 
there exists a constant c^ > and an exponent cr{k) > such that 



(6.1) 

and 

(6.2) 



\Utrix, t)\\Ll^^ < 



\\Urel{x,t)\\ 



r2 



< 



CkilogvY^'^^ 


yke 


Ckiiogvy^'^^ 



,ke 



\UtriO,t)\ < 



|Mref(0,t)| < 



Cfc(logt;)'^ 


k) 


yke 




CkilogvY 


ik) 



,k£ 



both uniformly on the time interval [^2,^3]. 

Let us first give an outline of the argument in this section. We would like to control 
||'U^||l°° l2, where w = u — Utr — Wrcf- If, after subdividing the interval [^2,^3] into 
unit-sized intervals, we could argue that w(t) at most doubles (or more accurately, 
multiplies by a fixed constant independent of q and v) over each interval, then we 
could take ^3 ~ logf and conclude that the size of ||w(t)||2,2 would only increment by 
at most a small positive power of v over [^25^3]- This was the strategy employed in 
[3]. The equation for w induced by the equations for u, Mtr, and Mref took the form 

(6.3) = idtw + |(9> - q6oW + F 

where F involved product terms of the form (omitting complex conjugates) w^, w'^Utj-, 
w'^Urei, wu^j., wutrUrei, and wu^^^, ti^ef'^tr, u^j.Urei, qSoUtr, qSoUrei- The Strichartz esti- 
mates were applied to (6.3) to deduce a bound on ||ty||iP ir for all admissible pairs 
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{p,r) over unit-sized time intervals [ta,tb]. Although a bound for {p,r) = (oo,2) 
would have sufficed, the analysis forced the use of the full range of admissible pairs 
[p, r) since such norms necessarily arose on the right-hand side of the estimates. 

A direct implementation of this strategy does not work for g < 0. The difficulty 
stems from the fact that the initial data w{ta) for the time interval [ta, tb] has a nonzero 
projection onto the eigenstate |g|^/^e~'''"^'. While the perturbed linear flow of this 
component is adequately controlled in Lr^ t ]^x^ ^^ ^^ equal to a positive power of |g| 
when evaluated in other Strichartz norms. For example, ||e2**'' |g|^/^e~'''"^'||2,6 le = 

\q\^^^. Thus each iterate over a unit-sized time interval [ta, tb] will result in a multiple 
of 1^1"^^^, and we cannot carry out more than two iterations. 

Our remedy is to separate from the above w = u — Utr — Urcf at t = ta the projection 
onto the eigenstate |g|^/^e~'^ll^'l, and evolve this piece by the NLSg flow. Specifically, 
we set 

Mbd(t) = NLS,(t) [P{u{ta) - U,,{ta) - U,,i{ta))] 

and then model u{t) as 

U{t) = Utrit) + Mrcf(t) + Mbd(t) + w{t) 

This equation redefines w{t) from that discussed above, and it now has the property 
that w{ta) is orthogonal to the eigenstate |g|^/^e~''''"^'. We will, over the interval 
[tajtft], estimate w{t) in the full family of Strichartz norms but will only put the 
norms L^ ^ iL^. and L'n' ^ ^H]. on Mbd(t). These norms are controlled by nonlinear 
information: the L^ conservation and energy conservation of the NLSg flow. The use 
here of nonlinear information is the key new ingredient; perturbative linear estimates 
for itbd are too weak to complete the argument. 

The equation for w{t) induced by the equations for n(t), Ubd(t), Utr(t), and u^dit) 
takes the form 

= idtw + \d1.w — q6o{x)w + F 
where F contains terms of the following types (ignoring complex conjugates): 

• (delta terms) Wtr^O; Uj-cfSo 

• (cubic in w) w^ 

• (quadratic in w) w'^Utr, w^Ureu w'^u^d 

• (linear in w) Wufj., WU^^^, WU^^, WWrefUtr, WUhdUtv, WMbdWrcf 

• (interaction) WbdW^r^ MbdM^f ^rdul^, UreiU^^, Wtr^bd' ^*trM?ef 

The integral equation form of w is 

wit) = Ugit)[{l - P){u{ta) - Utrita) - U,,i{ta))] 

- I f U^'^'-''^''' PF{t') + Uq{t - t')(l - P)F{t')\ dt' 
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We estimate w in the full family of Strichartz norms, and encounter the most adverse 
powers of |g| in the PF component of the Duhamel term. Of the terms making up 
F, the most difficult are the "interaction" terms listed above that involve at least one 
Ubd- Let A = \\u{ta) — Utr:{ta) — ^iref(^a) ||l2 • Then we find from the L^ conservation 
and energy conservation (see Lemma 6.1) of the NLSg fiow, that 

Using these bounds, we are able to control the interaction terms Wtr^bd and WrefUbd as 

(see Lemma 6.3). The estimate of the interaction terms comes with a factor |g|^^^, 
and thus the bound on the interaction terms with at least one copy of Ubd is of size 
A + \q\2A^. We want this to be at most comparable to A, so we need |g|2A^ < 
1. But the estimates of Phase 2 leave us starting Phase 3 with an error of size 
|g|3t;"6(^"^). Let us assume, as a bootstrap assumption, that we are able to maintain 
a control of size |g|3t;~6(^~2e) qj^ ^j-^g error. Then, even in the worst case in which 
A ~ \q\3v~6^^~'^'^\ the condition |g|5A^ ^ 1 is implied by f > |g|T4(-^+^'^). So, if we 
impose the assumption v > |g|T4(^+^^), we can carry out the iterations, with the error 
at most doubling over each iterate. If it begins at |g|3t;~6(^~'^), then after ~ elogf 
iterations, it is no more than \q\3v~6^^~'^'^\ the size of the bootstrap assumption. 

Now, to prove the error bound of size |g| 3t)~6(^~^) in the Phase 2 analysis required 
the introduction of the cubic correction refinement; the shorter argument in [3] would 
only have provided a bound of size |g|6t)~2+^. The bound \q\(iv~2+^ combined with 
the condition |g|2v4^ ^ 1 gives the requirement v > |q'|~^^, which is unacceptable 
since the most interesting phenomena occur for |g| ~ f . This is the reason we needed 
the cubic correction refinement in Phase 2. 

We now give the full details of the argument. We begin with an energy estimate 
for our differential equation: 

Lemma 6.1. If u satisfies 

idtu H — d^u — qSoU + \u\'^u = 0, 
then 

(6.4) \\dM^,t)\\Ll<'2\\d,u(x,0)\\Li + 2\q\\\uix,0)\\Li + \\u(x,0)\\l2. 

Proof. Multiplying by dtu, intergrating in space, and taking the real part, we see that 

Re I - / diudtudx — qu{0,t)dtu{0,t) + / \u\'^udtudx 1 = 0. 
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Here we multiply by 4 and integrate by parts in the first term: 

-dt I \d^u\'^dx-2qdt\u{0,t)\^ + dt j \u\'^dx = Q 

Integrating from to t and solving for ||92;M(t)||^2, we find that 

\\dMxMll= ||9,n(x,0)||i2+2g|n(0,0)|2-2g|n(0,t)|2 

+ \\u{xMU-\HxM\'Lt- 
Dropping the terms with a favorable sign from the right hand side, we see that 

WdM^Mli < \\dMx,0)\\l2+2\q\\u{0,t)\^+\\u{x,t)\\l, 

Next, using ||m||^4 < ||M||L2||M||i,g°, together with ||ii|||oo < ||u||i2 ||9^n||^2, we have 

\\dxu{x,t)\\l2 < \\d^u{x,0)\\l2+2\q\\\u{x,t)\\L2\\d:,u{x,t)\\Li 
+ \Hx,t)\\l2j\d^u{x,t)\\L2 
Here we use ||M(x,t)||L2 = ||m(x, 0)11^2: 

\\dMx,t)\\l2 < \\dMx,0)\\l2 + (2|g|||n(x, 0)11^2 + ||n(x, 0)||i2) ||9,n(x, t)|U2 

Using ab < ^a^ + ^b"^ to solve for ||(92;m(x, t)||, we obtain 

||5xM(x,t)||i2<2||9,M(x, 0)11^2 + (2|g|||M(x,0)|U2 + ||M(x,0)||i2)'. 

This implies the desired result. D 

We now give an approximation lemma analogous to that in Phase 1, but with the 
error divided into two parts. 

Lemma 6.2. Suppose ta < tb and tb — ta < ci. Let Ubd(t) be the flow o/NLSg with 
initial condition 

(6.5) Mbd(^a) = P[u{ta) - Utr{ta) " Mrcf(^a)], 

and put 

(6.6) W = U-Utr - Mrcf " "bd • 

Suppose \q\^\\u\^d(ta)\\L^ < 1 and, for some k eN, 



w{ta)\\Ll+c{k){qf^^^ + ||nbd(ta)|Ui + |g|(g)^hbd(ta)||i2 < C2 



\2 

V 

Then 



klU- , ,n < C3 [ \\w{ta)\\Ll + c{k){q) 



+ ||Mbd(4)||Li + kl(?)^ll«bd(ta)||i2 

The constants ci, C2, and 03 are independent of q, v, and e. 
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Once again we apply our lemma before proving it. Suppose now that for some large 
c and /c, where c is absolute and k depends on e, we have Iq'Kq')^ ||ty(t2)P < c and 
c(fc)(g)^ ^°^J.e — < c||t(;(t2)||L2. These additional assumptions cause the conclusion of 
the lemma to become ||w||l°° l2 < c||-u;(ta)||L2, and they will follow from assuming 
|g|io^q,^3^-i4(i-2e) ^ ^ ^^^ (g)t)~" < Ce^n, where c is an absolute constant and Ce,„ is a 
small constant, dependent only on e and n, which goes to zero when e — > or when 
n —>■ cxD. Let m be the integer such that mci < elogf < (m + l)ci. We apply Lemma 
6.2 successively on the intervals [^2, t2 + ci], . . . , [t2 + {m — l)ci, ^2 + mci] as follows. 
On [^2, ^2 + Ci], we obtain 

Applying Lemma 6.2 on [t2 + Ci,t2 + 2ci] and combining with the above estimate, 

ll'U^(',^)ll/'°° r? < Co||w(-, t2)llL2- 

II \ ' '"^[t2 + l,t2+2]-^x — Jll V ' ^n\J^x 

Continuing up to the fc-th step and then collecting all of the above estimates. 



tn,t 



2.'^3J 



Notice that Mbd(^) and w{t) are redefined by (6.5) and (6.6) as we move from one 
interval Ij = [^2 + (j — l)ci, ^2 +jci] to the next interval /j+i = [^2 + jci, ^2 + (j + l)ci] 
in the above iteration argument. That is, on Ij, Mbd(^) is the NLSg flow of an initial 
condition at t = ^2 + (j — l)ci, and on Ij+i, Mbd(^) is the NLSg flow of an initial 
condition at t = ^2 + jci that does not necessarily match the value of the previous 
flow at that point. In other words, at the interface of these two intervals, 

lim Ubdit) ^ lim u^dit) 

t/t2+jc^ t\t2+jci 

It remains to prove Lemma 6.2. On the way we will need to estimate the overlap 
of Ubd with Utr and Urcf. 

Lemma 6.3. Let ta < h, let Utr and Uj-ef satisfy (6.1) and (6.2), and let Wbd be the 
flow under NLSg of an initial condition proportional to e'^'^L Then 



(6.7) . _i 11 

<c |g| ^\\ubd{ta)\\LZ + \q\^{q)^\\uhd{ta) 



\h 



Proof. We prove the inequality only for 1 1 Utr^ibd 1 1 l°° l2 , the argument for 1 1 Urcf^ibd 1 1 l°° l^ 
being identical. We will use the following formula for Wbd^ 



ft 

„2, ,1 



Ubd(x,t) = c„e^*'' |g|2e''l^l+ / e-'^'-''>^^\ubd{x,s)\\M{x,s)ds, 
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where c„ is the overlap of Ubd with the hnear eigenstate, and is a constant bounded 
by the L^ norm of Mbd- We then have 

'l^trMbdlUff ,^,L2 < c„|g|2 ||nti.e^l''l||^^ ^^ 



[ta.ttl^ 



"[ta.tt] 



t2 



-i{t-s)H, 



''|Mbd(s)pMbd(s)rfs 



r oo r 2 



For the first term it is enough that the hnear eigenstate has L} norm proportional 
to |q'|~2. We remark that a better estimate is possible using the fact that Wbd and Wtr 
are concentrated in different parts of the real line. 



Ut 



,(x,t)|V'?l^'ldx < c\\utr{x,t)\\l^\q\-^ < c\q\-^ 



Putting back in the factor of ||Mbd(^a)||L2 |g|2, we obtain the first term of the desired 
estimate. 

Now we treat the integral term. 

r oo r 2 



Utr I e 

t2 

< c 



''^'-'^"^\uUs)\'uUs)ds 



t2 



T OO J 2 



where we have used the explicit formula for Wtr to take its supremum. We now use a 
Strichartz estimate: 



Here we use 



< c 



I'^bdl r4/3 ,1 + P|Mbd| rl r2 



< 



— \H\ II bdiii.[t^_t^]i.^ 
<c|g|^||nbd||ijy^^^^j^|||nbd||Lf^^,,^jLso. 

\\f\\Ll\\df\\Li. 

< c|g|2||Mbd||2oo i2||5xMbd|lioo ^2 

Applying (6.4) gives 

1 5 / 1 1 1 3 

< C|g|2||libd||loo r2 \\dxUhd{x,ta)\\l2 + H'AWhdix , t a^Li + ||tibd(a;, ta) 111,2 

But the L'i norm of Mbd is conserved, so that ||Mbd||L°° l"^ = ||wbd(^a)||L2- Moreover, 
because Mbd(^a) is proportional to e'''^', we have |9a;Mbd(a^, ^a)| = |?||wbd(a;, ^a)|, giving 
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||<9a;tibd(2^5^a)||L2 = ^l ||^bd(^a) ||l2 • This allows US to conclude 
i - / i - 1 - 

= C|g|2||^bd(ta)|ll2 [\q\^\\Uhd{ta)\\l2 + \q\mUbd{x,ta)\\l2 + \\Ubd{x,ta) 

^X \ ^X ^x 

Combining terms, we obtain 

<c|g|^(g)5||nbd(ta)||i2, 
giving us the second term of the conclusion. 

Now we proceed to the proof of Lemma 6.2. 
Proof. We begin by computing the differential equation solved by w. 

= idtu H — d'^u — q5o{x)u + u\u\'^ 
= idtw + -dlw - q6o{x)utr - q6o{x)urcf 

+ {Utr + Mrcf + Mbd + w) jwtr + Mrcf + Mbd + w|^ 
- I^tr^^itr - l^ircf^^ircf- l^bd^^ibd- 

This can also be written as the following integral equation: 



3 



u 



(6.8) 



w{t) 



-i{t-t2)Hq 



W(t2 



-i{t-s)H, [g5p(a;)^^^(5) + g5o(x)M,ef 



t2 



|Mtr(s)|Vcf(s) 



\w{s)\'^w{s)\ds. 



w||l- , ,Li + \\w\\Li , ,L-, and accordingly 



For this estimate we will use \\w\\x 
estimate the Ln ^ ,L^. norms of each of the terms on the right hand side of (6.8) for 
(p, r) equal to either (oo, 2) or (4, oo). We choose these norms in order to be able to 
apply the Strichartz estimate. 

First we treat the term arising from the initial condition. Here we use the fact that 
by construction, our initial condition satisfies Pw2{ta) = 0. 



\ r' T 



. < c\\w{ta)\\Ll 



Second we treat the delta terms, for which we obtain the following: 



<^-^)H^q5,{xW{s)ds 



<c(g)2||Mtr(0,t) 



tV Tr 






e-(*--)^,g5o(x)nref(s)rfs 



-^L.tt]-^^ 



<c{qf\\u,,f{0,t)\\L^ 

[ta 



[ta.t^] 



4/3 



We have used here the bound on tb — ta to pass from the LJ ^ , norm to the L|^ 
norm, and we have combined all the terms under the one with the least favorable 



tb] 



FAST SOLITON SCATTERING BY ATTRACTIVE DELTA IMPURITIES 



35 



power of {q). These last expressions are estimated using (6.1) and (6.2) respectively, 
to obtain, in both cases 



<c{q) 
d 

^''\w{s)\''^w{s)ds 



ke 



V 

Third we the treat term which is cubic in w. For this we have 

^-i{t-s)Ha\„..t „\\2 



T P T r 



< C I IIW^ILe/s ^<i/5 + \\Pw'^\\li Lr 






[ta.ti,] 



< C I l|u;||^8/5 ,18/5 + Ikilia r3r 



|5/9|| ||4/9 



At this stage we use ||w|Li8/5 < ||ty||/2 ||w;||/oo, and an analogous inequality for ||u;||i3r. 
We then use the boundedness of tf, — ta to pass to the L|^ ^ .Ll and Li ^ ,L^ norms, 
giving us 



< c\\w\ 



X- 



Fourth we treat terms of the form w^Wbd? as usual ignoring complex conjugates. 
This time we use p = 1, f = 2 on the right hand side of the Strichartz estimate. 

JP Tr 






We have used ||P/||l'-i < c|g|''2 '■i ||/||l'-2 to pass to the L^ norm, incurring a penalty 
no worse than Ig]^^^. This allows us to put an L^ norm on u^d, which is our preferred 
norm for this part of the error. 



< c|g|2|kbd|U- , ,L?lkl 



X- 



Here we use kh ||nbd||L°° l^ < 1: 

<c\Mx- 

Fifth we treat terms of the form w'^Utr and w'^u^ci- These terms obey the same 
estimate, and we write out the computation in one case only. 



<i--W^nj{s)W{s)ds 



tP Tr 






< c llii; Ut 
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The first term is the same as in the Wbd case. For the second term we pass to the L^ 
norm, so as not to be penalized in |g|. 



< C (||w||^||Mtr||L- , .Ll + IkVrlLi , ,L- ) 



Here we use the fact that our explicit formula for Utr allows us to control all of its 
mixed norms. 



< c\\w\ 



X- 



Sixth we treat terms of the form wul^, ^""ref; ^^"^ wux^t-u^^u again writing out the 
computation in one case only. 



<^~^)H.w{s)ul{s)ds 



T P J r 



< c I \\wu^ 



\ Va,th\ Va,th\ I'a.ibl ^ 1*0,*;,] ^ / 

< C(ta - Uf'^ {\\w\\l?? ^ .Ll + ||w||l4 r«> ) ||Mtjicx> rj^ 

<c{U-taf'^\\w\\x. 



Seventh we treat terms of the form wMti-Wbd and wu^^iUToA- 



-i{t-s)H, 



''w{s)Utris)Ubd{s)ds 



ta 



fP Tr 






\\PwUtrUbd\\a ,,Ll 



< C(l + |gh)||w||ii r^llUtrlU- .L-ll^bdlU- ^ .L^ 

3 1 

< C(t6 -ta)*|g|2||nbd||L-' , ,L2||w||x- 

< C{tb-ta)^\w\\x- 



where in the last step we again use \qU ||Mbd||L°° l2 < 1. 
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Eighth we treat terms of the form wu^^. 



<'-^)"'>w{s)ul^{s)ds 



^\ta.t^]^"^ 



<C\(1\\\WUI^\\ 4/3 1 

<c\q\\\w\\4/s ^Jlul^h?^ LI 

< c{U-ta)^2\\w\\x 



Ninth we treat terms of the form u^j.Urc{ and Mtr^ref) once again writing out the 
computation in one case only. 



-'^'-'^''^uUsfu.Mds 



^\ta,%]^"^ 



< C [ \\ulUrei\\r-i/S .i + ||PMt,Mrcf ||l1 , .^ 

< c\q\\\ulUrc{\\L'r- , ,Ll 
<c|g|||MtrMrcf||L-_,^jL2||Mtr||L-_,^jLi. 



At this point we use the exphcit formula for Ut,- to bound ||utr||L°° l2, and we split 
up ||MtrWrcf||L°° L^ into rcgious aloug the positive and negative real axis. 



< eg Mtr LfT , ,L^ Frcf Lff , ,L2^„ 



Ftr Lff , ,L2^J Frcf LP? , ,Lg° 



, (log^)-W 
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Tenth we treat terms of the form u^^ui^, ^^bd^tr^ref and Ubd^ref- 



<'~^^''^uUs)u,,{sfds 



,L'' 



< C (|kbd<|Li^,,^jL| + ||P«bd«?jLf,^,,L, 



<c|g|2||MbdMtrllLf» ,,L2 
< c|g|5||MbdMtr||L-^^^jL2 



II? 



< C M|Mbd(4)||Li + kl(9)^ll^bd(^a) 

Eleventh, and last, we treat terms of the form w^d'^tr and vt^^u^d. 



~'^'~'^''^uUs?uMds 



^\ta,H]^"^ 



< C ||MbdWtr||i4/3 1 + \\Pu^^Utr\\n ^ ,L; 

< c\q\\\ul^Ut,\\L^ ri 

< clglllMbdlUff , ,L2||MbdMtr||Lf» , ,L2- 

But using once again |q'|2 ||nbd||L°° ^2 < 1 we see that this obeys the same bound as 
the previous term. 

Combining these eleven estimates, we find that 



wWx < c \\w{ta)\\Li +c{k){q) 



2(logt;)'^W 



+ ||^bd(ta)||L2 + \q\{q)^\\Uhd{ta)\\l2 



+ ih - ta)^\\w\\x + \\wfx + Ikllx)- 



Choosing ci sufficiently small allows us to absorb the linear term into the left hand 
side. A continuity argument just like that in Phase 1 now gives us the conclusion. D 

7. Phase 2 for g > 

In this section, we prove Theorem 3. We repeat the three lemmas of Phase 2 in the 
case q > 0. These results are used in the course of Phase 2 in the case g < in order 
to approximate the flow of e~^^^° by that of NLSq, and they also give an improvement 
of the asymptotic result in [3] for g > 0. 



1/2 



Lemma 7.1. Let (p E L^ and < th. If q >0 and t^; < Ci\\(j)\\J, th 



en 



(7.1) ||NLSg(t)(/.-e 

in the Strichartz estimates. 



itHa 






<C2tf||0||i., 



where P and R satisfy p + -^ = |, and ci and C2 depend only on the constant appearing 
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Proof. Let h{t) = NLSq{t)(j), so that 

idth + -dlh - q5Q{x)h + \h\^h = 0. 

We use again the notation ||/i||x' = ||^||l°° l^ + ||/i||^6 ^e. We apply the Strichartz 
estimate, which in this case reads 

II "h>-tbr^ - ^ii^ii^ Ho.t^]^-^ 

once with (p, r) = (oo,2) and (p, r) = (6/5,6/5), and once with {p,r) = (6,6) and 
{p, f) = (6/5, 6/5). As before, the h terms are each estimated using Holder's inequality 
as 



II'^II^IS/S ,18/5 ^ ^t '""" 

yielding 



,18/5 ,18/5 < Ct^ ||«-||X') 



x'<cm\L^+C\\hrx> 

Using the continuity of ||/i||x'(t5)! we conclude that 

U'<2c||0|U2, 



so long as Sc^tJ/^ II 011 2 2 < 1. 

We now apply the Strichartz estimate to u(t) = h(t) — e"**^''^, observing that the 
initial condition is zero and the effective forcing term —\h\'^h, to get 



\\h{t) - e-^*^'0|Lp .« < C\\\h\'h\\^e^s ^e/5 < Ct'j'\\h\\\, < Ctf ||0||i.. 

Lemma 7.2. Under the same hypotheses as the previous lemma, 

(7.2) ||NLS,(t)0 - ^IUf^,,^,Li < c (t^ll^lli^ + tf 110111. + t,||0||i.) , 

where 

g(t) = e-''^''(f)+ f e^^(*~')^^|e-*'^^0|V^'^^0ds. 
Jo 

Proof. A direct calculation shows 

h{t) - g{t) = [ e-*(*~^)^^ {\h{s)\^h{s) - |e-'^^'0pe-'^^'0) ds. 
Jo 

The Strichartz estimate gives us in this case 



D 



S-ffo J, 1 2 „—isH„ 



\h-g\\L^^^Li < |||/i|^/i- |e-^^"'0|^e-*^"'0||^a/s ^6/5. 
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We introduce the notation w{t) = h{t) — e **^«0, and use this to rewrite our 
difference of cubes: 

\h\^h~\e 



-isHn ±\'2 „~isH„ 



^''^'>S\w\'' + e''^'Sw^ 



-''"^(t)\^w+ {e-''^^ ''^ 



w. 



-- w\w\ + 2e' 
+ 2|e" 

We proceed term by term, using Holder estimates similar to the ones in the previous 
lemma and in Phase 1. Our goal is to obtain Strichartz norms of w, so that we can 
apply Lemma 7.1. 

We have, for the cubic term. 



1/2 I 



W II 6/5 6/5 < CV \\w\\l?° , ,L| 



\w\ 



[O.tb] 



'^[ta,*^]'^- 



<ctt 



9 

L2- 



For the first inequality we used Holder, and for the second Lemma 5.1. Next we 
treat the quadratic terms using the same strategy (observe that as before we ignore 
complex conjugates): 



|e"^'*-^''(/)|wPL6/5 .6/5 < Cty^llwllioc r2||w||i6 L6 II 6 '"'■^"^ 



,6/5 ^6/5 



"lta,tb] ^' 



^fta.tl,]^^ 



< ct 



3/2 



7 
L2- 



And finally the linear terms: 



|||e-^'0p^||,6V5^^^6/5 < ctf |klUjf^,,^jL2||e--^'0||i6^^^^^,6 < ct,||0||i.. 
Putting all this together, we see that 



\h — qllroo t2 < clt 



imi^+C 



L2 +h 



D 



Finally we estimate, at time t2, the error incurred by dropping the integral term: 
Lemma 7.3. For ti < ^2 and = u{x,ti), we have 

(7.3) /%-^(*2-^)^^|e-*'-^'0|2e-^'-^«0rfs < 



Li 



{t2-ti) + {h-ti 



qe 



+ q^e 



2„-2v^ 



3 — 3i>l~'' 



where c is independent of the parameters of the problem. 

Proof. We write 0(a;) = 4>i{x) + 4>2{x), where 0i(a;) = e~**^*' /^e**^/^e*^''sech(x — Xq — 
vti), and estimate individually the eight resulting terms. We know that for large v, 
02 is exponentially small in L^ norm from Lemma 3.1 of [3]. This makes the term 
which is cubic in 0i the largest, and we treat this one first. 

I. We claim /^*' e-'^*''-'^"''\e-''"^(j)i\'^e-''"''(pids < c{t2 - ti). 
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We begin with a direct computation 

-t2 



Jti 



Ll 



< {t2-ti) ||e-*(*2-')^'|e-^"-^'0i|2e-*'^'( 



"^[?i,t2l^- 



isHa 1 \2 ^—isHn 



<c(t2-ti) |e-*^^'0ire 



'[tiM 



Li- 



lt remains to show that this last norm is bounded by a constant. This follows exactly 
the same argument as that given in part I of Lemma 5.3, with the difference that 
terms involving P are omitted. 

II. For the other terms we use simpler Strichartz estimates. The smallness will 
come more from the smallness of (/)2 than from the brevity of the time interval. 

-12 






Ll 



<c|||e~*'^'0i|V^'-^^ 



'2'lHi..l^i- 



We have used the Strichartz estimate with (p, r) = (1, 2), and use Holder's inequality 
so as to put ourselves in a position to reapply the Strichartz estimate. 



<c(t2-ti)5||e- 



1*1. '2I 



-isHo 



1*1. '2] ^ 



<c{h-h)"^\\MhU2\\Ll 

Note that here 6 = 1 — e. Similarly we find that 

ft2 



Jti 



< c(t2-ti)2g e 



2-2v^ 



Ll 



t2 



-*(t2-s)/f,. -is//, , |2g-i.H, , ^^ 



:^„3„-3?;^ 



< c{t2-ti)2q'^e 



Ll 



D 



-{1-e) 



These three lemmas improve the error in the case g > from v 2 to f 
in the language of [3], from w~2 to v~^. 
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